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Using the Rossi X-ray Timing Explorer we have measured lags of the 
9 to 33 keV photons relative to the 2 to 9 keV photons in the timing 
noise between 0.01 and 100 Hz in the accreting neutron stars 4U 0614+09 
^ ; and 4U 1705-44. We performed similar measurements on the accreting 

black hole candidates CygX-1 and GX 339-4 as a comparison. During 
£5 ! the observations these sources were all in low (hard) states. We find 

phase lags of between 0.03 and 0.2 radians in all these sources, with 
a variation in frequency much less than expected for a lag constant in 
time. We also measure a coherence consistent with unity in all sources. 
As already noted for the black hole candidates, these data are inconsis- 
tent with simple Comptonization models invoking a constant time delay. 
Comptonization in a non-uniform medium can perhaps explain the lags. 
However, the magnitudes of the lags imply that the hot electron gas 
extends to more than 10 3 Schwarzschild radii. This may constitute an 
energy problem. We argue that while a large hot cloud is possible for 
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black holes which may hide some of their accretion energy in advection, 
such a distribution may not be possible neutron stars where all the ac- 
cretion energy is eventually released at the neutron star surface. This 
casts doubt on the Comptonization model, though the energy problem 
may be resolved for example by a wind from the inner disk. 

Subject headings: accretion, accretion disks — black holes — stars: indi- 
vidual (Cyg X-l, GX 339-4, 4U 1705-44, 4U 0614+09) — stars: neutron 
- X-rays: stars 
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1. Introduction 

The X-ray properties of black hole X-ray 
binaries in the "low state" are strikingly sim- 
ilar to those of atoll sources (accreting low 
magnetic field neutron stars) in the "island 
state" (van der Klis 1994a,b). Both emit X- 
ray spectra that are dominated by a hard 
power law component, N(E) oc E~ a , whose 
photon index a increases (i.e. whose spec- 
tra soften) as the mass accretion rate, M, 
increases (Tanaka & Lewin 1995; Barret & 
Vedrenne 1994, van Paradijs & van der Klis 
1994). The power density spectra of the X- 
ray flux variations of both types of sources 
are dominated by a broad-band component, 
which is flat below a break frequency, V\>, and 
turns to a power law, P(y) oc v^ 13 above u^, 
with (3 ~ 1 just above the break (see van der 
Klis 1995). 

As M increases, so does z/ b ; the observed 
values of v\> are similar for the black hole and 
neutron star systems and cover the range 0.01 
to 30 Hz. As z/b varies, the high-frequency 
part of the power density spectrum tends to 
remain the same (Belloni & Hasinger 1990; 
Miyamoto et al. 1992), although this rela- 
tion is not always exact in detail. Generally 
it is observed that the power level of the flat 
part of the power density spectrum is anti- 
correlated with (and M) in black hole can- 
didates (Mendez & van der Klis 1997 and ref- 
erences therein) and in neutron stars (Yoshida 
et al. 1983, Prins & van der Klis 1997, Ford 
& van der Klis 1998, Mendez et al. 1998). 

Another characteristic of the low states 
of black hole X-ray binaries (BHXB) is the 
phase lag of high-energy photons with re- 
spect to low-energy photons ("hard lags") in 
flux variations covering time scales between 
~ 10~ 2 to ~ 10 2 s (see Miyamoto et al. 



1992). This phase lag is roughly constant 
over this range, at a level of order 0.1 radi- 
ans, and does not vary much with v\> (i.e. M) 
for a given source or between sources. These 
phase lags have been interpreted (eg. Payne 
1980) as travel time differences of photons in 
a hot Compton scattering electron gas (e.g. 
Sunyaev & Triimper 1979), which is also the 
medium in which the power law photon spec- 
trum is thought to be formed (e.g. Hua & 
Titarchuk 1995; Dove et al. 1997). In accret- 
ing neutron stars, lags have been measured in 
the kilohertz oscillation signals (Vaughan et 
al. 1997; 1998). These very small lags put 
stringent constraints on the physical size of 
any cloud which might Comptonize the pho- 
tons in the fast signals. 

It is clearly of interest to know whether the 
similarities between island-state atoll sources 
and low-state BHXBs extends to their time 
lag properties. We present here the results 
of a time lag analysis of two atoll sources at 
frequencies below 100 Hz and compare the re- 
sults to an identical analysis of two BHXBs. 

In Section 2 we describe the observations, 
and data analysis. The results of our analy- 
sis are presented in Section 3. They are dis- 
cussed in the framework of the Comptoniza- 
tion model in Section 4. We conclude that 
this model has fundamental problems. 

2. Observations and Analysis 

We have used data from the Rossi X- 
ray Timing Explorer (RXTE) proportional 
counter array (PCA) obtained during the first 
three observing cycles (see Table [I]). We 
choose four sources: CygX-1 and GX339- 
4, which are thought to contain black holes 
(Tanaka & Lewin 1995), and 4U 0614+09 and 
4U 1705-44, which are atoll sources showing 
X-ray bursts and therefore contain neutron 
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stars. In these observations the BHXBs are 
in the 'low' state and the atoll sources are in 
the 'island' state. They have similar Fourier 
power spectra which can be fitted by a bro- 
ken power law, with break frequencies of 0.18, 
0.07, 1.8, and 3.6 Hz for CygX-1, GX339- 
4, 4U 1705-44, and 4U 0614+09 respectively, 
and rms fractions of 29%, 36%, 22%, and 31% 
(2 to 9 keV, 0.01 to 100 Hz). The mean count 
rates are 4973, 586, 580, and 228 counts s" 1 
(2 to 9 keV, background subtracted for the 
whole PCA), with small long term variabil- 
ity. 

To quantify the difference in variability be- 
tween two energy bands we use the cross spec- 
trum defined as C(j) = Xl(j)*AT 2 (j), where 
X are the measured complex Fourier coeffi- 
cients for the two energy bands at a frequency 
Uj (van der Klis et al. 1987; for information 
on cross correlation analysis see Vaughan et 
al. 1994; Vaughan & Nowak 1997; Nowak 
et al. 1998). The phase lag, 0, between the 
signals in the two bands is given by the ar- 
gument of C (its position angle in the com- 
plex plane). The corresponding time lag at 
Fourier frequency, u, is simply (p/2ni>. We 
calculate an average cross-vector, C, by av- 
eraging over multiple spectra and binning in 
frequency, and then find 0. We calculate the 
error in from the observed variance of the 
C values in the real and imaginary directions. 
The resulting error is similar to, but slightly 
larger than, that derived from counting statis- 
tics (Vaughan et al. 1994; Cui et al. 1997), 
or (equivalently) from the coherence function 
uncorrected for counting statistics (Nowak et 
al. 1998). 

We calculate cross-spectra from data inter- 
vals of 256 sec and use a Nyquist frequency 
of 2048 Hz, which is high enough to avoid 
binning effects that dominate at frequencies 



above half the Nyquist frequency (Crary et 
al. 1998). As the resulting depends some- 
what on the choice of energy band, becoming 
larger for wider energy separations, we choose 
the energy bands to be closely similar. We 
choose two energy bands defined as effectively 
2 to 9 keV and 9 to 33 keV. These bands cor- 
respond to PCA channels 0-35 and 36-127 for 
the CygX-1 observations, and channels 0-25 
and 26-87 for the other data which were taken 
after a detector gain change in March 1996. 
We correct for deadtime effects by subtract- 
ing from C an average value obtained from 
800 to 1024 Hz where Poisson noise dominates 
(van der Klis et al. 1987). This correction is 
close to negligible. A positive value of cor- 
responds to a lag of the hard band relative 
to the soft band, a result we have confirmed 
by analyzing test signals and looking at data 
from the accreting millisecond X-ray pulsar 
SAXJ1808.4-3658. 

3. Results 

The results of the phase delay analysis are 
summarized in Figure [L|. In all four sources 
the hard photons lag the soft photons sig- 
nificantly. Even in 4U 0614+09, where the 
statistics are the worst, the lag is detectable 
with a significance greater than 5<r in the 1 
to 10 Hz range. Averaging over the 0.01 to 
100 Hz range, the values are 0.077 ± 0.003, 
0.093 ± 0.003, 0.092 ± 0.011, 0.067 ± 0.016 
radians for CygX-1, GX 339-4, 4U 1705-44, 
and 4U 0614+09, respectively. In 4U 1705-44 
and 4U 0614+09, the lags are significant only 
above 1 Hz; over the range 0.01-1 Hz is 
0.030 ± 0.014 and 0.047 ± 0.021 radians re- 
spectively. 

The significance of the phase lags is a fac- 
tor of 3 to 10 better for the BHXBs than for 
the atoll sources. In CygX-1 this is because 
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of the high count rate. In GX 339-4 it is a 
result of the high rms fraction of the noise 
and long integration time. We calculate the 
expected significances which are a function of 
the rate, rms fraction, observing time and co- 
herence (Nowak et al. 1998). Compared to 
4U 0614+09, the significances should be bet- 
ter by factors of 13, 14, and 3 in CygX-1, 
GX 339-4, and 4U 1705-44, respectively, sim- 
ilar to what we have observed. There is struc- 
ture visible for 6 function of frequency in 
the BHXBs (Miyamoto et al. 1992), in par- 
ticular increases somewhat with frequency. 
The lags are not, however, consistent with a 
constant time lag as shown by the dotted lines 
in Figure [TJ. 

In addition to the phase delay, we measure 
the coherence, defined by -y 2 = \ < C(j) > 
| 2 /(< IXxO')! 2 >< |X 2 (j)| 2 >) (Vaughan 
& Nowak 1997). The angle brackets repre- 
sent averages over multiple cross-spectra and 
Fourier frequencies. If the cross vectors have 
the same phase angles, 7 2 is unity; if the 
phases are random, r y 2 is zero. We corrected 
7 2 for the contribution from Poisson noise (see 
Vaughan & Nowak 1997); this correction sig- 
nificantly changes the calculated values. We 
calculate the error of the noise-corrected 7 2 
according to the prescription of Vaughan & 
Nowak (1997) in the high-coherence, high- 
variability regime applicable to our data be- 
low 100 Hz. Figure |2| shows the results. 
The coherence is consistent with unity for all 
sources. 

4. Discussion 

We have presented the first measurements 
of phase lags and coherence in the band lim- 
ited noise below 100 Hz in accreting neutron 
stars. We find that the phase lags of the atoll 
sources 4U 1705-44 and 4U 0614+09 in the is- 



land state are very similar to those of BHXBs 
in the low state as reported here and by pre- 
vious authors (Page et al. 1981; Miyamoto 
et al. 1988; Miyamoto et al. 1992; Nowak 
et al. 1998; Pottschmidt et al. 1998). In 
the noise at frequencies between 0.01 and 100 
Hz, hard X-rays lag soft X-rays. The aver- 
age phase delays that we measure in this fre- 
quency range are 0.09 and 0.07 radians for 
4U 1705-44 and 4U 0614+09, and 0.08 and 
0.09 radians for the low-state BHXBs CygX- 
1 and GX 339-4. The lags are not consistent 
with a constant time delay. 

We discuss our results within the frame- 
work of the idea that the X-ray photons in 
the power-law spectral component are the re- 
sult of Compton upscattering of low-energy 
seed photons, probably produced in the ac- 
cretion disk, by a very hot electron gas. This 
implies that the energy emitted in X-rays first 
resides in energetic electrons, which subse- 
quently lose part of their energy to low-energy 
photons in the scattering processes. In this 
Comptonization model, a low energy photon 
on average gains energy in each scattering 
event. Higher energy photons, on average, 
are the result of more scatterings. Therefore 
the higher energy photons emerge later than 
the lower energy photons produced simulta- 
neously, since the total path length they tra- 
verse is larger due to more scatterings. The 
measured unity coherence in all sources ar- 
gues for such a model as opposed to scenar- 
ios where signals at different energies are pro- 
duced in disconnected regions (Vaughan & 
Nowak 1997). 

In simple versions of Compton scattering 
in a uniform cloud, one expects a constant 
time lag between photons in different energy 
bands independent of frequency (Wijers et 
al. 1987), i.e. the phase lag, 0, is expected 
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to increase proportionally to the Fourier fre- 
quency, u, at least when the Fourier period 
is larger than the typical time delay This 
is inconsistent with the data as noted previ- 
ously. The recent work of Kazanas, Hua & 
Titarchuk (1997; KHT97) and Hua, Kazanas 
& Titarchuk (1997; HKT97), however, shows 
that a scattering medium with a radial den- 
sity gradient can produce phase lags roughly 
constant in frequency. 

We also note that in the atoll sources there 
are also quasi-periodic oscillation signals in 
excess of 1000 Hz which are thought to origi- 
nate very close to the neutron star. The upper 
limits to any hard phase lags in these signals 
is small (Vaughan et al. 1998). For a medium 
which induces constant time delays, the mea- 
sured lags of these fast signals should have 
been large. 

The size of the Comptonizing region is 
large. Our measurements of a 0.02 sec lag 
at 1 Hz for the neutron stars and 0.2 sec at 
0.1 Hz for the black holes correspond to a size 
of order 10 4 km for moderate optical depths. 
Even in the KHT97 and HKT97 models with 
non-constant density, the scattering medium 
extends up to at least 10 4 R S , where R s is 
the Schwarzschild radius. At the same time 
a large fraction of the total luminosity is re- 
leased in hard X-rays which are apparently 
produced by the Comptonizing medium. Rel- 
ative to the flux in the 0.1-100 keV band, in 
our observations Cyg X-l emits 65% of its en- 
ergy above 9 keV; in 4U 1705-44 this fraction 
is 25%. The question then arises: how can a 
substantial fraction of the emitted X-ray lu- 
minosity, which must originate from the con- 
version of gravitational potential energy into 
heat close to the compact object, reside in hot 
electron gas at distances of order 10 3 Rs away 
from the compact object? Energy production 



around accreting neutron stars was consid- 
ered early on (Zel'dovich & Shakura 1969). 
The energy problem discussed above has been 
raised by recent authors (e.g. Stollman et al. 

1987) . 

For the black hole systems one might argue 
that the energy budget in hot electron gas at 
distances of order 10 9 cm is sufficient to ac- 
count for a substantial part of the observed 
X-ray luminosity via Compton upscattering. 
The required mass accretion rate then has to 
be of order 10~ 7 M Q , i.e., several orders of 
magnitude above the value inferred from the 
observed X-ray luminosities of ~ 10 37 erg/s. 
The consequence of this is that the accretion 
flow within ~ 10 9 cm of the black hole must 
have a very low radiative efficiency, which 
suggests that it is advection-dominated (see, 
e.g., Narayan & Yi 1994). In this picture it 
would be the advection dominated accretion 
flow (ADAF) itself which provides the site for 
Comptonization. (Note that at 10 9 cm the 
virial temperature is already ~ 200 keV, in 
principle sufficient for the formation of a hard 
power law). 

However, a similar ADAF flow picture can- 
not apply to the atoll sources. These sources 
cannot have such a high accretion rate, since 
they have no way of hiding the correspond- 
ing photon luminosity liberated near the neu- 
tron star surface. The energetics is a prob- 
lem for disk configurations where the losses 
from Comptonization are expected to exceed 
the energy locally available from gravitational 
release or radiative heating (Shibazaki et al. 

1988) . KHT97 suggest that the Comptoniz- 
ing region is a quasi-spherical flow preheated 
by radiation from a central region (references 
in KHT97, Zel'dovich & Shakura 1969), per- 
haps a hot boundary region (Titarchuk, Lapidus 
& Muslimov 1998). Transporting a large frac- 
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tion of energy to large radii may also be ac- 
complished by a wind blown from the inner 
disk (Achterberg 1998). Magnetic fields of- 
fer another possibility for transporting energy 
from near the compact star to larger radial 
distances (Stone et al. 1996). To maintain the 
Comptonizing region in the neutron star sys- 
tems a large fraction (at least 25%) of the to- 
tal accretion power must be efficiently trans- 
ported from its point of release, within a few 
tens of km of the neutron star, out to dis- 
tances of order 10 9 cm. 

To summarize, we have shown that the 
phase delays and coherences measured in neu- 
tron stars are the same in X-ray binaries 
which likely contain black holes. If these 
lags are from Comptonization, the size of the 
scattering medium is large which may pose 
a problem for the energetics. Delays gener- 
ated by processes other than Comptonization 
may have to be considered as has been done 
recently (Bottcher & Liang 1998). The cor- 
rect model must describe both accretion onto 
black holes and neutron stars. 
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Table 1 
RXTE Observations 



Source 


Start Time 


Duration 


T 

- 1 res 


Nchan 




(UTC) 


(ksec) 


(//sec) 




Cyg X-l 


12 Feb 1996 12:45 


3.1 


16 


16 


GX 339-4 


3 Feb 1997 15:56 


13.1 


125 


128 




10 Feb 1997 15:40 


11.8 


125 


128 




17 Feb 1997 18:17 


7.9 


125 


128 


4U 1705-44 


1 Apr 1997 13:26 


12.0 


125 


64 


4U 0614+09 


22 Apr 1996 19:19 


26.9 


125 


64 




28 Aug 1998 19:30 


13.6 


125 


64 



Note. — Listed are the start time (in Universal Time, Coordinated) of each observation and 
the duration of data used. We use 'Event mode' data from the PCA with time resolution T res , 
and a number of channels N C han- 
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Fig. 1. — Phase delays between signals in the 2-9 keV and 9-33 keV bands for two black hole 
candidate binaries in the low state (left) and two neutron star binaries in the 'island state' (right). 
The dotted lines show constant time delays of 10 msec. Measurements which are less than la 
significant are plotted as 95% confidence upper limits. 
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Fig. 2. — Poisson noise corrected coherence functions between the energy bands 2-9 keV and 9-33 
keV. 
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